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Abstract: In this paper, the theoretical expressions of the specific heat and thermal conductivity of wood are 

first derived from the micro-structure of wood by applying some basic principles in statistical thermodynamics and 
physical mechanics. Tested by a great many experiments, the average error of theoretical values is about 5%, 
much less than that of all kinds of exiting empiric equations in the world. 
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Introduction 

Thermal properties of wood such as specific heat, 
thermal conductivity and so on are parameters which 
is required in scientific research and heat treatment 
of wood. Because wood is a kind of natural polymer 
organism whose chemical element and microstruc¬ 
ture are very complicated, it's very difficult to study its 
thermal property theoretically, therefore, for almost 
one century, it was by means of experiment that the 
scientists at home and abroad studied the thermal 
property of wood directly (Dunlap 1912; Kollmann 
1968; Cheng 1985; Gao 1985; Zhang 1986) But the 
question is that the empiric equations of the wood 
thermal property derived from experimental data by 
means of mathematical fitting have different func¬ 
tions. and different adaptable area (species tress, 
moisture content, temperature and so on), with a 
larger error, lacking the theoretical foundation. In or¬ 
der to solve the above-mentioned problems, the uni¬ 
versal theoretical expressions of wood thermal prop¬ 
erty are derived in this paper by means of basic 
principles of various natural subjects, such as statis¬ 
tical mechanics, physical mechanics, wood physics 
etc. 

Theoretical expression of wood specific heat 

The theoretical model to develop specific heat of 
wood is based on cell structure of wood, that is, the 
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pyran-rings, the basic unit of wood cellulose, are 
connected into a long-chain high polymer by glucose 
linkages. One pyran-ring has 20 atoms, 60 degrees 
of freedom, among which 54 degrees are of internal 
vibration including the stretching and bending vibra¬ 
tion of C-C, C-H, C-O, O-H bonds, and they all be¬ 
long to the Einstein’s mode of vibration. As to the 
other six, three of them are of translation freedom in 
the mass center of pyran-ring, the other three are of 
rotary freedom. These six degrees of freedom belong 
to the Debye’s mode of vibration. There are also 
oxygen atoms among rings, because the mass of 
oxygen atom between two rings is far below the 
whole mass of the atoms in the ring, it can’t form 
Debye's mode of vibration. So these three degrees 
of freedom belong to the Einstein's mode of vibration 
too. Therefore, the theoretical expression of specific 
heat of dry wood is given as: 
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acteristic temperature and Debye function respec¬ 
tively, and 0, is the characteristic temperature of 


each chemical bond among rings, R is universal 
constant of gases, p is molar mass of the glucose 
group. 

In order to simplify calculation, we can expand ex¬ 
pression (1) into Taylor series at T=273.15K, then, by 
substituting known quantity, we can obtain 


Co(t) - 1.2294 + 0.006714 t (kJ -kg _1 -K' 1 ) (2) 


It is obvious that there is a simple relationship be¬ 
tween the specific heat of dried wood C 0 and that of 
the wet wood C w (whose moisture content is w%): 
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Where C' is the specific heat of water. 

We calculate the specific heat of 33 species of 
wood with different moisture contents and tempera¬ 
tures by means of equations (2) and (3). When the 
theoretical values are compared with the experimen¬ 
tal values in the same conditions, it has been dis¬ 
covered that the maximum error is 7.8%, and the 
average error is 2.5%. 


Theoretical expressions of thermal conduc¬ 
tivity of wood 


In order to derive the theoretical expression of 
thermal conductivity of wood, we shall assume that 
the phonon is the carrier of energy in heat conduc¬ 
tion of wood. The phonon shows the particle property 
of lattice wave, and move irregularly in all directions 
in wood at sound velocity of u as gas molecules do, 
colliding each other and exchanging their energy. A 
phonon with velocity u, total energy E and mean free 
path L can transport energy E when it passes 
through the distance of L. 

Considering the double direction movement of 
various phonons, we find that energy is transported 
from the region of larger E to the one of smaller E, 
that is, from the region with higher temperature to the 
region with lower temperature. We can thus derive 
the theoretical expression of thermal conductivity of 
dry wood, and find 
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Where M is molar mass of wood cell, p is mass den¬ 
sity of dry wood, and other sings mean the same as 
stated above. 

The mean free path of the phonon L is related 
not only to the direction of heat current, but also to 
the frequency of sound wave in wood, temperature 
and density of wood, etc. Therefore, it is very difficult 
to decide the expression of the mean free path L of 
the phonon theoretically. According to the experi¬ 
mental results and theoretical analysis. We assume 
that there exists the relationship between the mean 
free path L and temperature T density of wood, De¬ 
bye characteristic temperature ©„ as follows. 
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Where p A is density of air (values in normal state), 
(7 and C 2 are two constants to be decided. (Cj is of 
length’s dimension) and can be decided according to 


the experimental data of the thermal conductivity. On 
the basis of better accuracy, we can use (7=10.15 x 
10 ~ 10 m, (7=7.987 (for heat current in the choral di¬ 
rection), and (7=12.60x 1CT 10 m, (7=5.283 (for heat 
current in the radical direction). 

In order to simplify calculation, we can expand ex¬ 
pression (4) into Taylor’s series at 1=273.1 5K, and 

expand T^ in expression (5) into power series of t. 
Then, by substituting known quantity, we can obtain 
separately simplified formula of wood thermal con¬ 
ductivity in the choral and radial direction. 
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The relationship between thermal conductivity of 
dried wood and that of wet wood (whose moisture 
content is w%) is shown as the following equation 

X w = (l +1 ,25w%]X 0 (8) 

We have calculated the thermal conductivity of 23 
species of wood in choral direction and that of 18 
species of wood in radial direction. It is easy to see 
that the maximum error of the former is 11.9% and 
that of the latter 13.7%, and the average error of the 
former is 4.8% and that of latter 5.3%, as the theo¬ 
retical values is compared with the experimental 
values. 
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